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LINEAR METHOD FOR THE CALCULATION OF THREE-DIMENSIONAL, 
LAMINAR AND TURBULENT BOUNDARY LAYERS 

H., W. Stock 
Dornier, Inc. 


1 . Introduction 


/I* 


Methods of calculation for three-dimensional, laminar and 
turbulent boundary layers and their application for three- 
dimensional aircraft wings and bodies will be described in this 
study. The methods of calculation are based on integral 
methods which, when compared to difference methods, suffer a 
loss of accuracy and are less flexible with regard to the use 
of various turbulence models. On the other hand, they require 
much shorter calculation times and storage facilities. 

The integral method for the calculation of turbulent flows 
is a‘n improved version of an existing method, the method for 
laminar flows is a new development. 

2, General Fundamentals 

2.1 Choice of Coordinates 


There are two directions marked for three-dimensional 
boundary layer currents, as shown in Fig. 1. One direction re- 
sults from the projection of the flow line at the outer edge 
of the boundary layer onto the body surface and indicates the 
direction of the main .flow t; the other one runs vertical to 
the main flow direction t and represents the crossflow direc- 
tion n. The distribution of the three-dimensional velocity 


*Niniibers in the margin indicate pagination in the foreign text. 
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profile into main flow and cros-sflow directions is sHown in 
Fig. 1. 

Three-dimensional boundary layers in the flow' line coor- 
dinates t and n were often calculated; where t and n are curvi- 
linear, orthogonal coordinates. Significant drawbacks result 
from the use of flow line coordinates. 1. The course of the 
flow li-nes at the outer' edge of the boundary layer must be 
calculated. 2. The fTow line coordinates' must be newly deter- 
mined for the same body geometry when the flow situation 
changes ' (change of angle of attack for instance). The intro- 
duction of arbitrary, curvilinear nonorthogonal coord-inates x 
and y removes* these drawbacks. In Fig. 1 the distribut-ionrof 
the three-dimensional velocity profile into the x- and y-direc- 
■tions is shown. There are simple relations- between the' . 
velocity profiles in the t 7 , resp. n- and x- resp. y-direction, 
in Fig. 2, through the angles a and X: 


U Si'n CX^-cL) -^Vcos CX-<^) 

I U. - ; 

I Si'n X 


Sin CX-oL) 
Sin X 

UsindL -h V cos 

Sin X 



V, 



sin <L 
sin A 


1 


I 


( 1 ) 


2 


Dkksinsb pxge is 

QB POOR QUAIiTX 



Ug is in this case the resultant velocity at the outer 
edge of the boundary layer, ui and Vi are its components in 
the 2 C- and y- directions. U and V are the velocity components 
in the t- and n- directions, see Fig. 1. 

A is the angle between the x- and y- direction; a is the /5~ . 

angle between the flow line projection at the 'outer edge of , 
the boundary -layer onto, the body surface and. the x-’directio,n 
and a+3 is the angle between the wallflow line and the x- 
direction. 

2.2 The Metric Coefficients of the Coordinate System 


When the body surface, for which the boundary layer is to 
be calculated, is given in Cartesian coordinates X, Y', Z »!thenl^ 
a curvilinear coordinate system x,y, can be so chosen on the 
surface that each point x,y defines a single point X,Y,Z. The 
transformation of the Cartesian coordinate system into the 
curvilinear one can be written as 

X = X(x,y) 

Y = YCx,y) C2) 

Z = Z(x.,.y) 

and the metric coefficients hi,h 2 and g of the curvilinear 
system are given by 
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^ ■DX QX , 'DYDY . 

0 X ^ 

* 

From eq. (2) one gets, after differentiation over x and 7, /_4 


, ^ f)X Q"x + 7)V ?)^Y . 

.’dx Ox 

(4) 

and 


, 2^!^ ^ Ox a^'x '■, O'y •a^Y ^ 

^0^ *2) X 0x *0^ Ox 0xOy Ox. OxO>^ 

with similar expressions for and 

lO X Q> y '~^x .'?) ^ 

2.3 Basic Equations 

It is sufficient for the derivation of integral equations 
to start with the equations of the laminar boundary layers, 

V » 

since the integral equations for laminar and turbulent boundary 
layers are identical. 
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In nonorthogonal cuvilinear coordinates the continuity 
equation and the x- and y- pulse equations are 


continuity equation 


O 


(5) 


x-pulse equation 

<? f IL ^ Jii. 

^ dx 

^ + a 

Qx ^ 'du 


i* w + CL^a,^ i- cl^v^ -h ci^uv^[. - 


'Cl^ 




Oa 

7>Vk 

02 


C6) 
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y-pulse equation 


\ 0X 


, 0 V 


+ w 


0v 

2>2 


0X 




02 




V- i>2 


£>3 a = 


(7) 


T and X are the shearstresses in the x- and y- directions, 

X y 

u,v, and w are the velocity components in the x,y and z direc- 
tion, ^ is the density and p the static pressure. The func- 
tions a- and b- are functions of the metric coefficients, as 
per appendix A. After integrating equations (6) and (7) from 
z=0 to z=6 (where 6 is the boundary layer thickness) and using 
equation (5) one gets, for the 


x-pulse integral equation 




l K, Ue '3x 9 ^ h,/ J , 


0 G>^i 

■ 


■( 8 ) 


/6 
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I +„© V ± 5 . (±i + a,l +A’' /-’4~' 

I^b-h, (A-Su ■J’SqiA,/ ■'Ih.U^'Ti^ a, 


^ A. 


fi d. 

t 


+ a, ~ -f Q-t — 

^ ^ a 


i 


+ a. 




y-pulse integral equation 

^ ^ ^ rf2-Me^; ^ 'Da 


K 


@..p 


^ o 


a ^ 




Dxlh,/ V /l; 


4) 0, 




C©' fOziV^i 2-/i.) + 6j + /i‘ Jll'^ 

h, Ue'^u Q.'dJhJ •• '" 


0= "9: 


■ 

C9). 

^ a 






h^ Ue Dy 


+ fc. ^ l ^@„b. 


I 


M is the Mach number o£ the outer edge of the boundary 

G 

and C£y are the componen-ts of 'the 


layer in this case and c 
■wa 1-1 friction coefficient ’ c^ in. the x- and y- 'directions ; 
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The equations employed here, which are known as "moment 
of momentum" equations in the English language literature and 
have so far been used only for the calculation of two-d!imensional 
currents, are derived in the following manner. The x-pulse 
equation, which is multiplied by the velocity component h, , 
is integrated from z=0 to z=6 and by using the continuity 
equation multiplied by u^/2 one gets the 


X- "moment of momentum" equation 


''D (I ) 'i 'DUi 

X I q 0 X 

■ oKiGTirKC-ymB 19- 

QE EQOR, QUALITY- 


ij.'dxl hj ^47 ^ 


a 



6 



©' 7 ^ 


V. 




■4- a -h Oi, — I -h I (© - (R) ) I — — -t a. ~ I 

" C/e" C/e" ] ^ ^ Ue 'Bo * C/e / 

, A ' ^ A. ? W . 


-25^ 


CIO) 


In the same way one gets the y- "moment of momentum" equa- 
tion by integrating the y-pulse equation multiplied by^hhe 
velocity component v from z=0 to z=6 and using the continuity 
equation multiplied by v^/2. 


y- "moment of momentum" equation /_^ 

I h, tjx '“A A, C4 ■ax Oo 


'A,. 


A 


I ^ f ^Us /?] , 7/ ? a' 






7 Ut * u7 i ^ A, cl 


2^ +i>, "■•’ 


'e tJ>c 


C/e 


( 11 ) 
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in addition the ’’entrainment” equation is used for calculation 
of turbulent flows since it describes the change of mass flow 
in the boundary layer. 

"Entrainment equation” 


( 12 ) 



The displacement thickness 6^ of the three-dimensional ^ 

boundary layer, which is often required as result of a boundary 
layer calculation, can be calculated from 


Z) 

^ ^ j 4 . j 

fSe ?V, cT'l 

1 _ ^ 9 ] 

1 7, 


* 

1 h. 7 




' h. / 


(13) 


Equation (13) has been derived by Myring .[1] according to 
the concept of the equivalent sources by Lighthill [2] . Equa- 
tion (13) can be solved as soon as the functions Ai and A 2 ’ are 
known. It should be mentioned that equation (13) is not re- 
quired for determination of the boundary layer development; 
it is solved independently of the boundary layer equations and 
serves only for the determination of the displacement thick- 
ness 6^. 
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The integral functions of the boundary layer in the x,y,z 
•coordinate 'System, Fig. 1, which appear in equations C8)-(13'), 
are listed in appendix B. The integral functions of the x,y,z 
coordinate system can be expressed by the appropriate integral 
functions of the t,h,z coordinate system, with the help of 
equation (1),- as per appendix C. 

Three-dimensional laminar and turbulent boundary layers in. 
curvilinear nonorthog'onal coordinates can now ‘be calculated 
with the.-,above listed basic equations. Velocity profiles for 
the main flow and cross-flow directions (t,n,z coordinate system) 
and temperature profiles are required for that. If those are 
available, the integral functions in the x,y, z .coordinate 
system can be calculated .in a simple manner, as shown in appen- 
dices B and'C. 

5. An -Integral Method for the Calculation of Three -'dimensional' /lO 

Turbulent, Incompressible Boundary Layers 

Cooke and Hall [3] have shown in their study that for three- 
dimensional boundary layers the flow in the direction of the 
main flow corresponds to a two-dimensional boundary layer flow. 
Myring [1] has accepted this assumption in his method for the 
calculation -of three-dimensional, turbulent, incompressible 
boundary layers through nonorthogonal,, curvilinear coordinates. 

He used the pulse equations (9) and (10) and the entrainment . 
equation (13) for determination of the boundary layer develop- 
ment. The description of the entrainment coefficient F, equa- 
tion (13) , is based on the concepts by Head [4] . The differ- 
ences between the present method and that of Myring are as 
follows: 

1. Power profiles of a single parameter for description 
of the velocity distribution in the direction of the main flow 
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have' been replaced- by Coles' profiles of two parameters [S] . 

The wallfriction parameter describes the velocity distribution 
close to the wall and the pressure gradient parameter deter- 
mines velocity distribution in the outer part of the boundary 
layer.' 

2. The additional empirical information required about 

the wallfriction coefficient, when power profiles are used, 
drops out iri this case. ' ’ 

3, Instead of the entrainment method used by Myring, in 

which the entrainment coefficient is determined directly, a 
"lag- entrainment" method has been employed in the present 
method which considers the nonequilibrium effects of the 
boundary layer, ’ , ' ' - 

5.1 Empirical Statements 

When calculating incompressible, adiabatic boundary layers /II 
the need ’to provide .temperature profiles for the determination 
of .in.tBgral functions is eliminated. The static temperature 
does not change in the boundary layer along coordinate z. 

5.1.1 Velocity Profiles 

3. 1.1.1 Direction, of Main Flow 


The velocity profiles with two parameters by 
Coles [5] are used, as described by the wall principle and the 
wake principle. 



k 






a 




C14) 
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where .. k is the von Karman constant, 6 is the boundary layer 
thickness,' tt- the pressure gradient parameter, W the universal ’ 
wake function which ‘can be approximated (6) through 




and ^ the wallfriction parameter which stands in the following 
relation to the Wallfriction coefficient c^^ 


T 



(16) 


The velocity distribution can also be written in / 12 
a form such as is used here for determination of the integral 
functions, with the- values for z=6, as determined from equa- 
tion (14) , entered into equation (14) xvith U=U. 

& 


=ffT[2-w(|j] -en|j 


3. 1.1. 2 Direction of Crossflow 


As in Myring’s [1] study, the statement by -Mager 
[7] 'and Johnston [8] describing the crossflow profiles as funcj- 
tions of those in the main flow direction, are applied here. 


Mager [7] suggested the following relation 


y 

Me 



(18) 
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where A has the following relation tO' the wall fl'ow line angle 
as seen in Figs. 1 and 2, which is equal to the angle between 
the projection of the flow line at the outer edge of the boun- 
dary layer of the body surface and the wall flow line, 

. A = tan $ . (19) 

• \ , J 0 hns 1 0 n [ 8 ] ~ jpr oy ide s jbr _the~"b oun dary Ta y e r • r. '> 

<fl^ next_ jBo~ the wall tlTe eorrelation - ' ' 

. tan 3 ■ ■ (20) /15 

e e ■ 

and for the remaining part of the boundary layer ' 

if- = A Cl-if-) C21) 

■ ■ e ■ e ■ . 

«. •- .. •- Johnston has shown that i.tr.is sufficient to 

use equation (21) for the calculation of integral functions 
of the crossflow. 

The relation between A and 3 is given by 
Johnston as ^ 

tan 3 = A (0.1 [c^^ cos 3]"^^^ - 1) (22) 

with the aid of the profiles in the directions of main flow 
and crossflow the integral functions in the t,n,z coordinate 
system can be stated as functions of 3,7r,a and A, in appendix 
D, which, in that way become variable functions in the proposed 
method of calculation. 


1'2 
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Integral functions of the x,y,z coordinate sys- 
tem are listed as ftinctions of a,X,.(S'ir,a and A in appendix E. 

The necessary derivations of the integral func- 
tions from a , A , 6 , IT ,!l? and A can easily be determined. 


5.. 1.2 Wallfriction. Equation /14 

So far only .three equations are available for- the 
calculation of four dependent variables 5,ir,^.and A; two pulse 
equations and one entrainment equation. A differential equa- 
tion which can b-e derived from^ the expres-sion for the imllfric- 
tion parameter T, will' be used as the fourth, inserting z=0 
with U=U in equation C14) results in an implicit equation .for 



Equation (23) , differentiated for x, represents the 
fourth differential equation that must be solved simultaneously 
with the pulse equations and the entrainment equation 



Si 

2>x 


? w 

'dx 


^ ^ ± = - 
^ Aq- 
• 


Ue. 


(24i) 


3,1.3 Entraiiunent Coefficient F 


Myring [1] used in his study the relation between F 
and the form parameter Hx as suggested by Head [-4] for two- 
dimensional flows. This empirical relation is based on mea- 
surements in boundary layers at equilibrium. Another sugges- 
tion was made by Horton [9], which takes into consideration 
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the so-called "upstream history effects" on the coefficient F 
in boundary layers of nonequilibrium. 

The entrainment coefficient P is not given as. the 
function of the form parameter Hi , as by Head, but is calculated 
through a differential equation, which is solved simultaneously 
with the. abovementioned four differential equations. 


,'2f C2S) 

’ '^X ' 

I ^ 

and where Fq is the entrainment coefficient under conditions of 
equilibrium, 6^ is a characteristic boundary layer th-ickness and . 
C is a constant. 


Horton' [9] assmes that .entrainment is closely bound 
up with the average shears tress in the outer part of the 
boundary layer’ and showed' that measurements in the boundary 
layers at equilibrium can be correlated with 


' F 

< 'e 


0.1 zz 


(H^ - 23) 




(26) 


where ^ represents the shearstress for z/6‘=0. 5. A mixing 

method relation of the Prandltl type is used for conversion of 

* 

equation (^26) into a relation between and C8U/-8 zO‘q 5, 



where l/'X = 0.083. 


(273 
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With 'the aid of Coles’ profiles it can.' be shown that 

w 


V 


u. 


4 s? 


^ ' os 


H, - a*3 


(28) 


The equations'. (26) - (28) result in- 


with 



= ytSSS 



5e ^eV 




0.63 f 

f 


(29) 


(30) 


In Horton’s study [9] the characteristic boundary 
layer thickness 6^ in equation (25) was equal to the pulse los 
thickness 6ii and .the’ corresponding value for C was C=0.012. 
it turned, out, however, that boundary layer thickness 6 repre- 
sents a more suitable measure of the characteristic thickness 
since it fits- better into the order of magnitude of the big. 
"eddies" that produce entrainment. .The corresponding value 
for C is C=0.1. This C value corresponds to a typical relax- 
ation length of 10 and agrees with the value by Bradshaw [10] 
in the equation for scaling of the dissipation length, which 
takes on..a<iform similar to equation' (25) . 

Since entrainment occurs in the outer part of the 
bqundary layer, where velocities of the crossflow are small, 
the two-dimensional view can be extended to three-dimensional 
flows when entrainment is coupled to the velocity profiles in 
the main flow direction,. 



5.2 Final Form of the Differential Equations 
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Equations (8), (9), (12), (13) and (25) can be written 
with the aid o£ the empirical statements 

X’pulse equation 


0 M +< 5-0 ^ 


(31) 


y- pulse equation 



OcT 

9x 


+ 


5r 

'dx 


+ 

^Z'l 


JJ- 


+ 

Z'i 




= J, (32) 


1 


entrainment equation 


■0 ^ 
^3 "ax 


4- 




2,r X 


-f- 


<S0 


or 


•^ 7 T c^x 


cT^ 





(33) 


Wallfriction equation 



T^J 


■f 2c^ 





= i 1(34) 

I 

1 


entrainment coefficient equation 


9F 

ax 



(35) 
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displacement area equation 
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8 6 ^ 
9x 


= ]6 


(36) 


The expressions are listed in appendix F and the func- 
tions are listed in appendix E. 


5,5 Numerical Integration 

Equations (51) -(35) are solved simultaneously. If desired, 
the displacement area of the three-dimensional boundary layer 
can be calculated with equation (56) . 

Myring [1] has pointed out that the differential equations 
,(51) -(35) are of hyperbolic nature. Of the resultant three 
characteristic directions the two outer ones describe the areas 
of influence and dependence .in a calculated point. These areas 
are limited by two- strai'ght_lines that pass through the calcu- 
lated point and form the angles go" and a+g with the x-a-xis. 

When determining the y-derivations contained in the functions, 
this situation is taken into account. • 

When both angles a and a+3 are positive backward^/ directed 
differences are used 


y 


‘ Lf ' m.ti 


— Qrn,n Qm,n-- 


(37) 


when a and a+3 are negative forward directed differences are 
used 

- Q. 
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(58) • 
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and when a and a-*-3 have different signs centered differences 
are used 



Q 





y 


(39) 


Q is a dependent variable, m is the 'counter for integra- 
tion in the x- direction and n the counter for . integration in 
the y-direction. 


When the y-derivations of the dependent variables along 
line x=constant are known integration in the x-direction can be 
accomplished through an explicit intermediary step method 


Q 


Q 


m,n 


-J- AX 


0Q\ 

/ m+ d 


(40) 


^fx^m+1/2 n found by extrapolation of the value (Q) m+l/2,n 
at a distance of Ax/2 from (Q)^ .with (-|$)„ _ 

JIl j H O ^ Hi y 

^H^-H/2,11 calculated with CQD„^.i/ 2 ,n' 


Boundary conditions are required where characteristic di- 
rections run into the calculated area from the outside. For 
the three-dimensional boundary layer this corresponds to a 
situation in which the boundary layer flow enters the calcula- 
tion area from outside. It also follows that boundary condi- 
tions are not required in places where the boundary layer flow / 20 
leaves the calculation area. 


For the initial conditions of the calculated results shown 
in 3.4 the measured values of 9‘n, H and 3 are used and. 
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converted into the dependent variables and A, as per 

appendix G, The initial value for coefficient F was determined 
from the measured value of form parameter Hi and the relation 
valid for flows in equilibrium., by Head [4] . 

5.4 Comparison with Experiments 

5.4.1, van den Berg and Elsenaar [11] 

Berg and Elsenaar ’s experiment is of special signifi- 
cance since it best simulates the flow conditions at sweptback 
wings of all known measurements of three-dimensional boundary 
layers. Berg and Elsenaar have measured the turbulent boundary 
.layer, on a yawed plate (angle of sweepback 35f) in the low 
velocity area. The wall of the i^ind cha'nnel above the panel 
was so constructed that the generated pressure increase was 
large enough to cause separation of the boundary layer. Beyond 
that .it has been attempted to approach conditions for an infin- 
itely long sweptback’ wing, through shaping of the eiid plates of 
the panel measured. 

A Jionorthogonal, straight line coordinate system' was 
used' for the calculation in which the x-direc-tion is identical 
with the approach flow direction and the y- direction is paral- 
lel with the leading edge of the panel. 

In' Fig. 3 the results for the moment loss thickness 
01U3 the form parameter H of -the velocity profile in the main 
flow direction and the wall flow line angle 3 were .compared 
.with the measurements. In addition the results of calcula- 
tions, according to the methods of Cousteix [12] and Smith 
[12] , are shown. Cousteix used the pulse equations and the 
entrainment equation.- He determined the velocity profiles 
from the similar solutions of the turbulent boundary layer. 
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Smith expanded Myring:'s [1] method to compressible 
flo'ws and carried out the verifications shown here of experi- 
ments in the incompressible area. The results correspond, 
therefore, to Myring’s method. The curves by 'Smith, shown in 
Fig. 3., correspond to the results shown at the Euromech 60. 
The;-.leit'fe.a::s :M and J in Fig. 3,., and in all subsequent figures, 
next to the present method and that of Myring, correspond to- 
the results for crossflow profiles by Mager or Johnston.: For 

the calculations in Fig. 3 the mea's.ure'd wall pressure distri*- ’ 
bution and the condition of infinitely long, yawed wings was 
used as input. Up to a panel depth of x~0.9 the present pro- • 
cedure- shows better agreement with the measurement than that of 
Myring . 


"Small digressions in the distribution of the resul- 
tant velocity U and of angle a resulted from' different inter- 
pretations of the data measured , in Fig. 4. In. Fig, 5 the 
results of calculations for crossflow profiles' according to 
Mager are shown for input of different velocity data at the 
edge of the boundary layer. The curves '/!/ correspond to the 
measured values pf U and of a, the curves /2/ were determined 
from U- and the condition of infinitely long yawed wings and 
curves /'3/ from .the ^^all pressure distribution measured and 
the condition of i'nfintely long yawed wings. From the results 
It becomes clear that small disturbances in the velocity data 
of the friction-free flow for boundary layer flows close to 
detachment, can provide widely varying results. Fig. 6 shows 
the corresponding results, as in Fig. 5, for Johnston cross- 
flow profiles. In Fig., 6a a comparison with difference 
methods (Schneider [12], Krause [12]) is demonstrated. Curves 
/!/ and /3/ correspond to the conditions mentioned before. 

No advantage can be seen in the more complicated difference 
■method. The detachment of the boundary layer is not ' 
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calculated by any difference method. For crossflow profiles, 
according- <to Mager and Johnston, the pres.ent method deter- 
mines the point of separation in the .area of the experimen- 
tally determined separation, after input of measured values 

for and a', 
e 

5.4.2 Johnston [14] 


Johnston [14], examined the development of the boun- 
dary layer in the area of low velocity, as shown in Fig. 7. 

Air comes out of a rectangular channel and flows between two 
test plates, the height of the channel -corresponding to the- 
distance between them. The rear end of the measured distance 
is formed by a wall that is placed at right angle to the test 
panels and to the channel axis. The velocity distribution in' 
the plane that is. imagined stretched in the middle between the 
test panels, corresponds to the distribution for a two- 
dimensional ray incident on a wall at right angle. Since too 
few values are given for the outer velocities to calculate 
the boundary layer at the lower test panel, such velocities 
are .used as would result for such a configuration according 
to potential theory [15] . The calculation is carried out with 
orthogonal, line.ar coordinates, the results are shown in Figs. 

8. The calculated results by Smith [13] are shown for compar- 
ison. (A continuous curve of Smith's results, cannot be repro- 
duced since his results -were shown only for discrete points*.^ 

The differences between the results of the present method and 
that of Myring are only small. 

3.4.3 Vermeulen [16] /25 


Vermeulen [16] measured the development of the boun- 
dary layer at the bottom panel of a rectangular channel that 
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had a* 60° curvature-. Fig. 9 shows a sKe.tch o£ the measured 
distance together with the flow lines at the edge of the boun- 
dary layer and the- wall flow lines. A curvilinear coordinate 
system was used in which the x-direction followed the test 
stations in the flow di-rection, with the y-direction at right 
angle -to it. The values measured on lines A and E were used 
for limit .conditions. Figs. 10 show a comparison with measure- 
ments and- results by Smith* [13] . The digressions of the 
present method from the measurements are clearly smaller than 
those of Myring's method. In Figs. 11 and 12 the velocity 
profiles in main flow and crossflow directions, that were mea- 
sured and calculated for Mager crossflow profiles, are shown 
along line C. For the test- case of Vermeulen and for the tests 
(11), (14), Figs. 13-17, the calculated results from both 

•methods show the expected differences. In the area close to 
the wall and in the outer part of the boundary layer. Coles' 
profiles show better agreement with the measurements- for the 
main flow direction than the* power profiles of Myring. Agree- 
ment with the crossflow profiles is bad in both cases, though 
measurements are reproduced slightly better by the present 
method. 


• 5 . S Discussion 

The two parameter profile family by Coles, for description 
of the velocity profile in the direction of the main flow, 
showed better results for three-dimensional boundary layers 
than the- use of one parameter power profiles, something al- 
ready known for the two-dimensional case. An attempt has been 
made to include conditions of nohequilibrium, that should be 
considered for the calculation of turbulent boundary layers, 
by means of the "lag- entrainment”' method. Agreement between 
measurements and calculations is good considered- the effort 
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-invested in the method.- The proposed integral method provides 
good results even when compared with difference methods. 

4. An Integral Method for the Calculation of Three- 
dimensional, Adiahaticj Laminar , Compre'ssihle 
Boundary Layers 

The pulse equations (8) and (9) and the; "moment of iiiomen-^ 
turn"' equations (10) and (11) should b'e employed for- the calcu-- 
lation of three-dimensional laminar boundary layers. As in the 
case of turbulent boundary layers, chapter 3, temperature and 
velocity profiles are required for determination of the integral 
functions appearing in equations (8) -(11). 

4.1 Velocity Profiles 


It is assumed that velocity profiles for the main flow and 
crossflow directions can be determined from similar solutions 
of two-dimensional laminar boundary layers. Since the flows in 
question are compressible but adiabatic, 'a solution of the 
Cohen-Reshotko equations (17) can be dispensed with. The solu- 
tion of the Faulkner-Skan equation (18) for similar incompres- 
sible flows provides the foundation for the putting together 
of the profile families, together with the Stewartson (19)’ 
transformation which provides a correlation between the adi- 
abatic, compressible and incompressible flows. The use of the 
Stewartson transformation also alleviates the need to generate 
a- temperature profile. 

4.1.1 Similitude Solutions fOT the Two-dimensional 
Incompressible Boundary Layer 

A simple type of solution for the two-dimensional 
Prandtl boundary layer is the similitude . solution. These are 



.solutions that are so constituted that velocity profiles a-t 
various distances X- can be made to- fit them through proper 
scaling of the velocity U and the. distance at right angle to 
the wall,' Z (Z is identical with the independent variable y 
used otherwise for two-dimensional problems). For this case 
the boundary layer equations ar.e reduced from partial differ- 
ential equations to regular- ones . 

Similitude solutions exist for ’when the velocity of ‘ 
the potential flow U is proportional to a power of the running 
length measured from the stagnation -point. 

Ug(X) = const X™ • (41) 


The similitude transformation of the independent variable Z, 
which then leads to an ordinary differential equation, is 


n 


=-z( 


mi- 


-I Ue 




■ (42) 


The known Faulkner-Skan differential equation then becomes 


fi I . + £f i.t + 3„(;i-f >2) = 0 


(43) 


with limit conditions 


n = 0: f=f’=0 and’ n = oo: f» = l (44)- [1J_ 

df u 

where ^ ^ ^ xepres.ents the p,ressure gradient param- 

eter which stands^in the following relation to the exponent in 
equation (44) : 
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2m 

m ■+ 1 


C45) 


■ Fig. 18 demonstrates the results for U/U^ for 
2^3^- 0.199, The curve for 3=0 corresponds to the Blasius 
solution for a flat plate, 3= - 0.199 provides the velocity 
profile for retarded flow, which leads straight to -sepanafion, 
and 3=2 is valid for highly accelerated flows. 

4.1.2 Direction of Main Flow 


' It is assumed that the velocities shown in Fig. 18 
can represent the velocities in the direction of the main flow. 

4.1.5 Direction of Crossflow /28 

It is assumed that the velocity profiles in the direc- 
tion of the crossflow can be represented by 


V 

U 




e * 


( 46 ) 


V X 

where (ff-l describes a certain type of profile and c is an 
0 

arbitrary constant. A possibility for generatiqn of the var- 
ious types of profiles in the crossflow direction is shown in 
Fig. 19. The sectioned line represents the reference profile 
(U/U )^, which corresponds to the velocity, distribution U/U 
for 3=2.0 with a freely chosen value for n^6=6.75. n6 is the 
transformed boundary layer thickness and is defined as' the 
distance from the- wall in the transformed plane, for U-/U =0.99. 
The curves drawn correspond to velocity distributions U/U^ for 
values of 2.0 £ 3 _<-0.18, The following relation provides the 
curves for 
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^U. 


'U. 


( 47 > 


which includes also such profiles of the "cross-over" type 
V X ■ ' 

wher.e <rf-J changes signs in the boundary layer. T.o maintain 
-conditions of symmetry after the passage through "cross-over" 
profiles, curves III and IV correspond to the curves for $= 2.0 
■and 1.0 with changed sighs. Curves I and .II result from inter- 
polation between curve 11.1 and the curve' for B'= - 0.18. ' 


4.2 Determination of Integral. Functions 
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The example of the energy loss thickness a-t the mai-n flow 
profile 



is t'o show how the integral functions of the profiles in the' 
main flow and crossflow direction, see appendix C, can- be de- 
termined with the aid .of the velocity profiles in Figs. 18 and 
19 and the Stewartson transformation [1.9] . The Stewartson 
transformation 


cL2 d 


o 


(49) 


permit,s the recalculation of compressible integral functions 
into incompressible ones. and a^ are the density and sound 

velocity in a reference state of the gas,.) Equations (48) and 
(49) lead to 
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where 6^ is the transformed boundary layer thickness in the in- 
compre.ssible .plane. 

Velocity profiles are shown as £un.ction of* the similitude 
variable n in Figs. 18 and 19. With the aid of equation (42)', 
equation (5.0) can be redefined to 


(5h) 


with 


The moment loss thickness of the velocity profile in 
direction q£ the main flow, 6ii, is chosen as scaling function 
for the boundary layer. 

As per equation (51) 


with 





K 


•i>i^ 


do dj 
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Equations (513 and (.52) then show 


0111 ^ 0ii4^ (53) 

-^1 1 

* ^ ■ * 

All phy-sical integral functions (6iii, etc.) can 

therefore be expressed by-the physical scaling function 0ii and 
by the rela-tion- of the corresponding transformed integral func- 
tions (Kill, Ki, etc.). 

h i ^ 

The transformed integral functions, which depend only on 
the profiles in the' main flow direction,^ shown in Fig. 18, are 
demonstrated in Fig. 20 (see appendix H for definition of the / 51 
functions) ; plotted over the parameter ' * " 


a 



Kii 


(54) 


where 


£ ’ * 

w 


f !•) 

^ dri '^w 


a=0 corresponds to the detachment profile 8= -0,199 and a-0, 38954 
to the profile of the. strongly accelerated flow 3=2.0. ■ 

. The transformed integral functions,, which depend only on 

V ^ 

the profiles in the crossflow direction (^) , see Fig. 19, are 
demonstrated in Fig. 21 (for definition of® the functions see 
appendix H) plotted over the parameter 


b - K[-~ 



(55) 
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K 2 is the transformed displacement thickness of the crossflow 
V X 

profile (ij-) for a transformed boundary layer thickness n6=l. 

The transformed integral functions K 112 , K 221 , Ki 2 and L 12 
(see appendix H for definitions), which depend on the profiles 
in the directions of main flow and crossflow, are shown in 
Figs. 22-25 plotted over b with a as parameter.. For determina- 
tion of these mixed functions the transformed boundary layer 
thickness of the crossflow profile has been set equal to the 
main stream profile. 


The following approach is suggested for determination of 
the constant c in equation (46) , which permits calculation of 
the desired crossflow velocity profile V/U from (V/U The 

physical displacement thickness in the crossflow direction is 





With equation (49) we get 


(56) 




Consequently constant c is 



( 59 ) 


6 ^ ±_ 

^ 7<r 

All physical integral functions in the t,n coordinate sys- 
tem can then b.e.*d.eterminod in the following' manner . 

' Q'(t,n)' =• f C0 i.iJgCc)h(a;b) , . (60) 

where Q' is an integral function. All the required physical 
integral functions have been listed in appendix H. 6n,a,b 
and c’ thus become the -dependent variables in .the method of /55 

calculation suggested .here . 

The physical integral functions in the x,y coordinate sys- 
tern, which appear in the pulse equations x and y used here- and 
in the "moment of momentum” equations (8) -'(^11), can be calcu- 
lated from the physical functions of the t,n. coordinate -system 
through the relations given , in appendix C. 

‘■Q(x,y) = Q(ei-i,a,b,c,.cf.,A-) ■ ' ' (61) 

The required derivations of functions Q after Oh i,.^a,b,c,a 
and’ X can easily be determined. . , - 



The following holds ' 



wherein Qq ' = - Ab ' ^ fojc instance. 

o 1 1 oti 1 1 
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The expressions are listed in appendix I. 

Equations (62) -(65) are solved simultaneously. If so de- 
sired^ the displacement area of the three-dimensional boundary 
layer can he determined with equation (66) . 

4.5 Numerical Integration 

Numerical integration is handled in the same way as ex- 
plained in section 3.3. The functions 6a i, H, and the wall- 
flow line angle $ are required as initial conditions. The 
transformed integral functions, which are given as functions of 
a or b, or a and b, are made available, through polynomial fits. 

For the start of the calculation at the stagnation point 
or stagnation line of a body, the sign of the wallflow line 
angle 3 decides the sign of the crossflow velocity (^) , Fig. 
19, for the first integration step. ® 
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At the start o£ the calculation only such profile types 
can be considered for which 


holds . 






I 



r 


Once the crossflow has developed from the stagnation pointy /'5S 
or stagnation line^then the profile type- of the crossflow 
changes through change of the sign of function i.e., the 

flow line at the outer edge of the boundary layer has a turning 
point. '’Cross-over'* profiles appear and after passage through 
this type of profile, profile forms as shown in III and IV, Fig. 

19,' are reached. Should a renewed change occur in the sign of 
the function , then "cross-over" profiles will reappear. 

The profiles will in that case pass through in the same direc- 
tion as during their first ' appearance . 


4.4 Comparison with the Calculated Results of a 
Difference Method 


The results of difference methods for laminar flows can be 
used for comparison in estimating the quality of approximation 
methods such as the integration procedure. The difference 
method was developed by Horton [20].. Three test cases were 
used for infinitely long yawed wings with 45*^ sweepback in the 
incompressible region. A rectangular coordinate system was 

IT 

used with the y-direction running parallel to the leading edge. 
The velocity component at the outer edge of the boundary flayer 
in y-direction, vi, is for all test cases. 



( 67 ) 
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with corresponding to the component o.f undisturbed initial 
flow velocity in the x-direction. .The velocity component* at 
the outer edge of the boundary layer in the x-direction, u-j , 


is 

Test case 1: ^ = 3(| - for | > 0 (68) 

^00 

Test case 2: ='l ' for | < 1.0 (69) 

- -- -- - = 1 - 0.567(^ - 1) for J > 1 (70) 

^OO ^ ^ 

Test case 3: - ^ = 1 for — < 1 (71) 

^ = 1 - 0.1134(| - 1) for I > 1 (72) 

OO 


c is a reference length-. Test case 1 corresponds to a boundary 
flow that starts at a stagnation line. The test cases 2 and 3 
are boundary layer flows which correspond at first to. a plate 
boundary layer and are then s.uddenly exposed to an Increase in 
pressure. The calculations were made for a Reynolds htimber 
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The comparison of results from calculations is shown in 
Figs. 26-28. All calculations end with the detachment of the 
‘laminar boundary layer. In test case 1 the sign of the wallflow 
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line angle changes but the displacement thickness of the cross- 
flow direction da remains positive, which means that there are 

■yr 

"cross-oyer" types of crossflow profiles in the range ^ for 
3>0. The present integral method gives a very satisfactory 
d'escription of these complicated flow conditions. Beyond that 
the separation point is determined accurately. For test case 
2 results are equally satisfactory, but small digressions occur 
in the determination of the separation point. During .recalcu- 
lation of test cases 1 and 2 agreement was looked for in the 
wallfriction coefficient c^ during establishment of the initial 
‘conditions and deviation in the form parameter 



was permitted. * For test case 3 two different results of calcu- 
lation are shown in Fig. 28, with agreement of C£ and H for the 
initial conditions. Here, too, the agreement is satisfactory. 

4.5 Discussion 


The results show that the idea for using the concep-t of 
the "moment of momentum" equations in the three-dimensional 
case 'as well, leads to good results. 

5 . Boundary Layers on Wings 

* 5.1 The ■ Influence of Wing .Tapering on the Development . . 
of a Three-dimensional, Turbulent Boundary Layer on 
a Transsonic Wing 

The development of a three-dimensional boundary layer on a 
wing was often calculated to an approximation with methods that 
were based oh the concept of the infinitely long yawed wing. 

In comparison to this quasi -two -dimensional way of looking at 
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•it one aspect o£ three -.dimensionality, the wing tapering, will 
be examined closer here. 


5.1.1 The Metric Coefficients 


To make the calculation simpler the boundary layer' . 
development . will be developed here not for the given wing con-, 
tour but for a flat* plate with' a top view corresponding to that 
of the wing on which the potentialitheoretical pressure distri- 
bution is' superposed . The calculation is carried out in a non- 
orthogonal linear coordinate . system. The x-direction runs 
parallel to -the direction of initial flow, the y-direction is 
identical with percentage lines on the. wing. 

Fig. 29 shows the ‘'G'a-n'ee.sdcan coordinates X,Y,Z and 
not the nonorthogonal, linear coordinates x,y for the infintely 
long yawed wing and .the tapered wing with ■ straight leading and 
trailing edges. 


The ■ following connection results 
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X = a(y) 

; ''' 

i Y. = y 


C73> 


2 = 0 


Th"e functions a(y) and bCy) .are shown- in. Fig. 29*. 

From equations- (3) and (4) we get the metric coefficients and 
their derivations to (also on Fig. 29) 


h. - 1 - - 


1 : 

^T< 
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- 0 
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; g = coi X Sin X^ ('I - 




(77D 


C78) 


<79) 


(80) ^ 


(81) 
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The angle A in a point x,y corresponds to the angle be-, 
tween the lines x=const. and y=const. at the point x,y. Aq is 
the angle between the coordinates x and y at the origin of 
the coordinates. " _ . ' • 


The following holds true 



.(8-3) 
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is the distance from the origin of the coordinates to the 
point where the extended leading and trailing 'edges of the wing 
touch. Cyx=‘» for the infinitely long yawed wing';) 

5.1.2 Description of the- Various Types of /41 

■ Calculations 

Fig. 30 shows the top view of the wing, the wing sec- 
tion examined and the sweepback of the leading and- trailing 
edges 'in that wing section. 

Fig. 31 shows the distribution of the resulting flow 

velocity for that wing at the outer edge of boundary layer U , 

dimensionless as the initial flow velocity over the wing- 

chord at the suction side in the section examinedt c is the 

wing chord of the wing section .examined. In addition, the run 

of angle is given. The leading edge sweepback of the wing in 

the section examined is (j)^=32°, the .trailing edge sweepback is 

(j)H=16®. The sectioned line for function a in Fig . - 31 .presents 
- " ' - » ' 

the run of a for an infinitely long yawed wing with a sweepback 
of /■2=24'’ for identical velocity distribution 

‘The same” integral method is used for two boundary layer 
calculations (13) . Case 1 corresponds to the. results for an 
infinitely long yawed wing with 24° sweepback. The wing for 
case 2 is tapered with a constant leading' edge taper of 32° and 
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a constant trailing edge taper of 16°, corresponding to the 
wing of the wing section examined. The distribution of and 
a over the. wing chord, hj percentage, is identical with that on 
the wing of the section ex-amined. 

5.1.3 Discussion of Results ' '/42 

The results for Mager- cross flow profiles are shown in 
Figs. 32-34. The moment "loss thickness 6ii increases slower 
in the trailing edge area, i.e., in the area of retarded ex- 
ternal flow at the -tapered wing, than at the infinitely long 
yawed wing. The same holds for the form parameter .H.. Since 
the form parameter H- increases with retarded external flow, 
the wing taper acts to weaken the retardation of the effective 
external flow. ... 

The wallflow line angle & is ,an effective measure for 
the three-dimensionality of the flow (a+g>0 means that the 
boundary layer material flows toward the wingtip'^-. Fig. 34 
shows clearly that tapering weakens the three-dimensionality 
of the flow. 

5.2 The Development of the Laminar; Boundary Layer on a 

Transsonic Wing 


The development of a laminated boundary layer on an. infin- 
itely loiig yawed wing in* the region of transsonic velocity is 
shown in Fig. 35. ■ The velocity data at the outer edge of ‘the 
boundary layer correspond to those that are also used for the 
turbulence calculation in section 5,1.2, Fig.. 31, Th.e calcula- 
tion is carried out with" a constant value c, equation (46) . The 
turbulent values of 0-iiiiand B are shown for cora,par ison in Fig. 
35. The moment loss thickness is. significantly smaller in the 
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laminar case. This result becomes understandable when comparing 
the development of the boundary layer on a flat plate in the 
laminar and turbulent case. For laminar flow 

011 ~ /43 


and for the turbulent case 


The velocity -decrease of the external flow for xs;0.5 is so great 
that laminar separation occurs. In the turbulent case no separ- 
ation ‘occurs here. It is known that turbulent boundary layers 
can tolerate greater pressure increases in the external flow 
before they separate. Results for the wallflow line angle 3 
show clearly that the three-dimensionality of the flow is far 
more pronounced for laminaT boundary layers than for turbulent 
ones. A similar result is also described in reference [22]. 

6. Turbulent Boundary Layers on Bodies 

[Revolution Ellipsoid) 

The development of a turbulent boundary layer on a revolu- 
tion ellipsoid (ratio- of main axes 4:1) has been calculated for 
an attack angle of 0® and 10®. 

6.1 The Metric Coefficients /44 


The connection between the Cartesian coordinates X,Y,Z 
and the coordinates of the’ revolution ellipsoid, as per 
Fig. 36 is 

X = a cos Tp 
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Y = b , s ill ’1’ 
Z=bsinycos^ 


C86) 


In this case a and b are the main axes of the revolution, ellip- 
soid. 


Using equations (3) and (4) for the metric coefficients 
and their derivation we continue 
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6.2 Results 


Tlie data for the friction- free flow around the revolution 
ellipsoid were determined through the procedure by K. Maruhn 
[22]'. Figs. 37 and 38 show the distribution of the resultant 
velocity of friction-free flow U and of angle ot for an attack ’ 
angle of a =1.0“', plotted over ^ and x/a. The distance x and 
the angle a are defined in Fig., 37. The angle ^ in the tangen- 
tial plane at a surface point P is defined as the angle between 
the vector U and the straight line that is generated by sec- 
tioning the tangential plane with a plane that stretches through 
the x-axis and the point P, as in Fig. 36. 

The. results of a boundary layer calculation for a Reynolds 
number 


Re = = 5.2 X 10® 

i^Ref 

< 

where the undisturbed initial flow is the reference condition, . / 46 
are shown in Figs. 39-41. The functions 0ii.,and>H are shown in 
comparison for an attack angle of a^-0“* 0.j i and with it all 
boundary layer thicknesses increase for the incident ellipsoid 
not only in the x-direction but also in the direction of the 
.circumference:. The. minimum* 'values lie in the plane of symmetry, 
on the side facing the wind. For value.s of-x/a>0.5 the wall- 
flow line angle on the topside assumes negative values while 
being positive on the side- facing the wind. In the plane of 
symmetry of top and bottom side there are no crossflows (6=0^ • 

The course of 0ii and particularly of -H for values of x/a>0.5 
is noteworthy. The maximum values of these functions do not 
lie in the plane of symmetry on the topside C^=180“) , but for 
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5P:2l50“. This may be explained throug-h the piling up of boun- 
dary layer material because o£ the cross-flow, which changes its 
sign at the circumference. G’eissler [23]- has observed- similar 
conditions in h.is study (development of the .laminar boundary 
layer for an incident ellipsoid) . The maximum increase of the 
form parameter H on the side ojf-the ellipsoid may also indicate 
the free separation of the -vorticity layer [24], which was 
examined closer by Wang [25] for an incident ellipsoid. 

7 . Summary 

Methods have been introduced for the calculation of three- 
dimensional, turbulent and compressible, iaminar flows. In the 
turbulent method the power profiles, which were used in Myring’s 
work for the description of the velocity profiles in main stream - 
direction, were replaced by Coles’ profiles.,. 'Be-yond that a 
"lag-entrainment" method was introduced instead of the simple 
"entrainment" concept. The quality of the results of the cal- 
culations., as compared to the experiments, could be improved 
even more if more suitable models were available for description 
of the- cross-velocity profiles. The proposed integral 'method 
gives good result's in comparison with difference methods. 

The laminar method,, which uses "moment of momentum" equa- 
tions for the first time for three-dimensional flows-, gives 
good results in comparison with an exact difference method. 

Calculations for wings showed that tapering of wings re- 
duces the three-dimensionality of the flow. For identical 
pressure distributions the .laminar boundary: layer separates 
earlier’ and shows a more pronounced three-dimensionality than 
turbulent flows. 
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The deyelopment o£ a; turbulent boundary layer was calcu- 
lated for an ellipsoid for 0° and 10 ° ' incidence . The boundary 
layer thioknesses for incident ellipsoids are smaller on the 
bottom side and greater on the topside than f.or ,an ellipsoid 
with 0° incidence. A noteworthy result is that the.-maximum 
values of boundary layer thicknesses and of the form parameter 
do not occur at the apex of the topside- but on the side oh the 
part of the revolution ellipsoid turned away from the wind. 
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Appendix A 
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Definition of functions a- and h. 
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Appendix B /55 . 

The integral functions in the x,y,z coordinate system 
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For laminar flows the dissipation integrals and can be 
written as 


I 

i 

1 




) 

I 


di \ 




o 


(B-2) 


]x isi’here-the dynamic \:iscosity of the gas 
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Appendix C 

Relation between the integral functions in the x,y,z and t,n,z 
coordinate^ system 


/55 


(Q = — 

Si 


+ Z S/ j Sln^CX-c^) cosCk-^) 
•h 3 022^ S/n CX-oL) cos^(X-<) - COS^X 


( 0 ) 




^ r ' 

wx[9.,s;n«c s.n^a-oc) + <j;« j s,oa-. 4;«sa-^) ' 
+ 0j2^s.Wcos‘-a-^) + e_^^cos<s,Vi'(:;L-^; ’ (c-d 
" ^ 0221 ^ g^coS-CCoS^CA-^jl 


0 


^2.-1 




a 


Sf'n ("A -C^) f + S/j Zs,'n^ CDS ^ 5/n CA 


■*■ ^22^ COS'^od s/n (X-fl*::) - sin V cos Cx~ 






-S/no^coscCCo^ ^TA-ccj - 0^^ oosV CO sCX~^jj 


&, 


2 ZZ 

K 


-d r t 

+ 2 SinUcos< -h 3 0 , 




' -.v^*,!. ;t 
Kk' . 

V 


Si^SS^^ page is. 

^ SUAiUCCa 


50 



Continuation: 
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Continuation: 
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^ " m ^ 

The integral functions appearing in equation (C-1), which are 
based on velocity profiles in the main flow and crossflow . di- 
rections, become ' - 
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Appendix D (Turbulent Boundary Layer ^ 
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Calculation of the integral functions of the boundary layer in 
the t,n,z coordinate system 

The boundary layer integral functions can be expressed as func- 
tions of the boundary layer thickness 6, of the pressure grad- 
ient parame-ter tt', the wallfriction parameter a and the parameter 
A, with the help of the -Coles' profiles for the main flow direc- 
tion and the statements of Mager or Johnston for the crossflow 
direction. We write f=f (it,0,A) . 
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For the integral functions which depend only on the velocity 
profile in the' direction of the main flow we write: 
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For the remaining integral functions: 
--Mager Crossflow profiles-- 
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Johnston Crossflow profiles-- 
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For the expression we get; 
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k=0.41 is here the von Karmon constant and Li and Lz calculate /62 . 
to : 



The derivations needed for the calculation 



can be determined from the above equations . 
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The integral functions o£ the boundary layer in the x,y,z. coor- 
dinate system. 

With the, aid of the relations between the integral functions in 
the t,n,z and x,y,z coordinate systems, as determined from appen- 
dix C and the f functions determined from appendix D, the 
integral functions in the x,y,z coordinate system can be named. 
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Continuation: 
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The f coefficients that appear in the equations (E-2.) are 
given in appendix -D. 


The derivations required for the calculations 



can be easily determined. 
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The right-hand sides of equations (31) -(36) 


The expressions I on the right-hand sides of the equations (31)- 
(36) produ ce _ 
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Continuation : 
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Appendix G 


(Turbulent Boundary Layer) 
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Calculation of the initial condition 


When recalculating the experiments the initial values of the de- 
pendent variables 6,ir,a and A, must be determined at the first 
test station from the measured values of Su, H and B. 


For the form parameter H we get 





and with equation ID -2.) 
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o; 
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In addition the following is true 


rf- = 1 H e.. 
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T ^ ^ ^ ^ A 

Equations (G--2) and (-G-3) are' iterated with an estimated value’ 
of IT. A can. be calculated for Mager profiles from 


A = tan B 


(G-4.) 
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and 'for Johnston profiles with equation (22) . 

The initial value for the entrainment coefficient F.-Kasuthe cor- 
relation by Head [4] . 


with t 

i 





! 


(G-5) 


IbRIGINAC PAGE IS 
DD io_Qg QUAIJX2 



Appendix H 


(Laminar Boundary Layer) 


/71 


The Transformed and Physical Integral Functions 
1. The Transformed Integral Functions 


Integral functions which depend on ‘profiles in the direction of 
the main flow (n^ is the transformed boundary layer thickness 
for which U/U^=0.99) 
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Integral functions, which depend on profiles in the direction of 
the crossflow (all functions standardized for rj6=l) 
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Integral functions -which depend on profiles in the main flow 
and crossflow directions (ri6 in the crossflow direction is set 
equal to the ri6 of the main flow direction) 



(H-11) 


(H-12J 
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. The Physical Integral Functions in the t,n, Coordinate 
System 

The definition of these integral functions is given in appendix 
C. Integral functions which depend on profiles in the main flow 
direction 
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CH-19) . 


Integral functions . which de.pend on profiles in the crossflow 
direc-tion 
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In-tegral functions which depend', on profiles in the main flow 
and crossflow directions 
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01 ij a, b and c are the dependent variables and 


1 Ri 


e. 


Ue Se 1 
/^e 


(H-31) 


i noe 


Me 


(H~32) 


I'^here t is the ratio of the specific heat of air. 
The Chapman constant c^ is defined as 
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and are ’the static temperatures of the gas at .the wall, 

resp. the outer edge of the boundary layer. 
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The right-hand sides o£ the equations C-62) - (66) 
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Fig. 1. Three-dimensional Boundary Layer Flow 
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Fig. 2. View of the Plane Tangential to the Body' Surface 
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Fig. 6, Resuits for Johnst6.n crossflow profiles fox input of 
various data of the potential flow 
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Fig. 10- Continuation 
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Fig. 13. Comparison of the measured and calculated velocity profiles in the maiii 
flow direction 







Fig. 15. Comparison of the measured and calculated velocity profiles in the 
' main flow direction 
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Fig.' 18, Velocity profiles (laminar^ in the main flow direction 
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Fig. 28. Comparison of calculation results of a difference 
method with the present integral method 
(Laminar boundary layer, test case 5) 
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Coordinate systems at tapered and untapered wings 













SIT 



Fig. 31. Distribution of velocity U^/Up^Qf and of the angle a over the wing chord 












Trace of the moment loss thickness over the wing chord 
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Fig. 34. Trace o£ the wallflow line angle over the wing chord 
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Distribution of the angle a at the incident revolution ellipsoid 
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Fig. 41. Trace o£ the wallflow line angle B for the incident 
revolution ellipsoid (turbulent boundary layer) 
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